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Abstract 

In this paper we consider thin films of AB block copolymer melts confined 
between two parallel plates. The plates are identical and may have a preference 
for one of the monomer types over the other. The system is characterized by 
four parameters: the Flory-Huggins parameter accounting for the repulsion 
between the A— monomers and the B— monomers, the fraction / of A— monomers 
in the block copolymer molecules, the film thickness d, and a parameter h quan- 
tifying the preference of the plates for the monomers of type A. In certain regions 
of parameter space, the film will be microphase separated. Various structures 
have been observed experimentally, each of them characterized by a certain sym- 
metry, orientation, and periodicity. We study the system theoretically using the 
weak segregation approximation to mean field theory. We restrict our analysis 
to the region of the parameter space where the interaction strength % is close 
to its spinodal value x,, the composition / is close to 0.5, the plate preference 
h is small, and the film thickness d is close to a small multiple of the natural 
periodicity Lq = In/qo (the first three conditions define the weak segregation 
regime). We will present our results in the form of phase diagrams in which 
A ~ I {d/nLo — l)/£ I is placed along the horizontal axis, and T ^ {% — X.v)/£^ is 
placed along the vertical axis, where £ is a measure for the composition asymme- 
try (/ — 0.5). We present a series of such phase diagrams for increasing values 
of the rescaled parameter H = he^^. We find that if A is small, corresponding 
to films whose thickness is commensurate with the bulk periodicity, parallel ori- 
entations of the structures are favoured over perpendicular orientations. We also 
predict that on increasing the value of the parameter H, the region of stability of 
the bcc phase shrinks. 
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1 Introduction 



Bulk systems of AB block copolymers may undergo a phase transition to a state which 
is homogeneous on a macroscopic scale, but which is phase separated on a meso- 
scopic scale. Many experimental 111131 and theoretical 0|ll|31IHI studies of bulk block 
copolymers have been published, enhancing our understanding of their phase behav- 
ior. More recently, there has been an increasing interest in thin block copolymer films. 
Such films can be created by spin-coating a block copolymer solution on a substrate, 
and then letting the solvent evaporate. One can create either a free film for which 
the upper surface is formed by air while the lower surface is formed by the substrate 
QIHIEIIIOI; or a confined film for which the upper surface is a solid plate parallel to 
the substrate [TT\. There is a qualitative difference between these two cases. If the film 
is confined, the microstructure might have to be compressed or stretched in order to 
fit inside the film. If the film is free, however, the system can avoid this distortion by 
the formation of terraces lfT2l l9l. so that locally its thickness is everywhere commen- 
surate with the natural periodicity of the structure. In this paper we will only consider 
confined films. 

Using X-ray scattering, various microstructures have been observed in copolymer 
films. The simplest of these is the lamellar structure EUTHIHl. The lamellae may 
be oriented either parallel, or perpendicular to the film. In lfT3l fT4l [TSl fT6l [TTl [TSl 
the relative stability of these two lamellar structures was investigated theoretically for 
the case of confined films. Also more complex structures have been reported experi- 
mentally 1 10 1, such as cylindrical structures, hexagonally perforated lamellae, and the 
gyroid structure. The cylinders of the cylindrical structure can be oriented either par- 
allel, or perpendicular to the film. The relative stability of these structures has been 
investigated theoretically in ifT^EHlUn . 

Apart from experimental and theoretical research, also some Monte Carlo studies 
have appeared ifTTllSm . In [[17 1 the phase behavior of symmetric AB-diblock copoly- 
mers confined in a thin film was investigated using both Monte Carlo simulations, 
and self-consistent field theory. The film boundaries were assumed to have a prefer- 
ence for the minority component. For thick films and small AB-incompatibility the 
authors found a lamellar structure in which the lamellae are oriented parallel to the 
film. For thin films and larger incompatibility, the authors find transitions between 
parallel and perpendicular orientations of the lamellae. In lEUI an A^B^gA^ triblock 
copolymer trapped between two plates having a preference for the end blocks was 
studied using Monte Carlo techniques. It was found that in very thin films, a cylin- 
drical structure with the cylinders perpendicular to the film arises, while in thicker 
films the perpendicular cylinders are replaced by either parallel cylinders, or parallel 
lamellae, of perforated lamellae. 

For experimentalists it would be useful to know in advance in which part of the 
multi-dimensional phase space one can expect to find interesting microstructures. This 
kind of information can be obtained using theory or simulation. We can distinguish 
two types of theories. On the one hand we have theories which require heavy numer- 
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ical calculations, like the self -consistent field theory. For a given point in phase space 
they will give an accurate prediction for the structure to be expected, but since these 
methods are rather time consuming, they are not suitable for obtaining an overview 
over the whole of phase space. On the other hand there exist the more simple (ana- 
lytical) theories like the strong segregation theory, and the weak segregation theory. 
Due to the approximations made, their predictions are less accurate, but they have the 
advantage that they can provide physical insight, and that the whole phase diagram can 
be mapped out rather quickly, providing a global overview. The rough results of such 
theories can then be used as a guide for more time consuming numerical methods or 
simulations, or as a guide for experimental work. 

2 Equivalence between a confined film and a bulk sys- 
tem 

We consider a thin film of molten AB—diblock copolymer trapped between two parallel 
imprenetrable plates. The copolymers are assumed to be flexible, in the sense that 
the persistence length is much smaller than the block length. Apart from the chain 
connectivity, the excluded volume of the monomers, and the excluded volume of the 
plates, which represent strong interactions, two weak interactions are present in the 
system. The first is a short-range interaction between the different monomer types, and 
is taken into account by means of the Flory-Huggins %— parameter. If the %— parameter 
is positive, there is a tendency for the A— monomers to separate from the 5— monomers. 
The second is a short-range interaction between the plates and the monomers. The 
range of this interaction may be either comparable with, or shorter, or longer than 
the size of a monomer unit, but it must be short compared to the radius of gyration 
of the blocks. If the interaction strength between the A— monomers and the plates is 
different from the interaction strength between the B— monomers and the plates, there 
will be an accumulation of one of the monomer types near the plates. The presence 
of these two weak interactions may lead to a non-homogeneous spatial distribution of 
the two monomer types. Our aim is to find the time-averaged shape of the equilibrium 
composition profile in terms of the parameters of the system. Having this goal in mind, 
it is convenient to switch from a microscopic description in terms of the positions of 
all individual monomers, to a mesoscopic description in terms of the coarse grained 
composition profile \|/ defined by 

^{x,y,z) = pA{x,y,z) - f (1) 

where pA is the coarse-grained A— monomer concentration, and / is the system 
averaged A— monomers concentration (we assume that the volume of a monomer of 
type A equals that of a monomer of type B). The units have been chosen such that the 
total density of the system equals unity, from which it follows that < / < 1. Note 
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that there are many microscopic states corresponding to a given coarse grained state 
\\r. We will analyze the system by calculating its mean-field free energy F. We make 
the assumption that this free energy can be obtained by minimizing a functional F[\\r] 
having the form 

= -5[X|/] + AF[X|/] (2) 

The functional F [\\f] can be considered as a generalization of the Landau free en- 
ergy known from the study of ferromagnetic systems. All effects concerning the weak 
interactions (monomer-monomer, monomer-plate) are present in the second contribu- 
tion AF, which vanishes if these interactions are absent. The first term 5[\|/] is the 
entropy of the coarse grained state \|/; that is, the logarithm of the sum of all a priori 
probabilities of the microscopic states which are compatible with this coarse grained 
state. It takes into account the effects of the strong interactions (excluded volume, 
chain connectivity), and it corresponds to the Landau free energy of a film in which 
the A— monomers are physically indistinguishable from the 5— monomers. It is pos- 
sible to write down an expression for S in the form of a series expansion in \\f. The 
coefficients of this expansion can be related to the single-chain correlation functions 
of the system in which the weak interactions are absent. We will assume that the 
third and higher order correlation functions can be built from the second order ones 
G(?',^i; 7,^2) (see [3 1). In principle, these correlation functions could be calculated by 
summing the a priori probabilities of all microscopic states in which, of a given chain, 
monomer / is at positional, and monomer j is at position ^2- Within mean-field theory 
it can be proven that the correlation functions thus obtained satisfy reflecting boundary 
conditions at the plates (see |22|, where this was established for the first time). For 
a half-infinite system occupying the part of space satisfying z > 0, the second order 
correlation function can be written as 

G(/,i?i; j,i?2) = Gbuik{ij Ri'Ji Ri) + Gbuik{h Ri'J: SRi) (3) 

where S represents a reflection in the plane z = 0. If equation Q is generalized 
to a film, the right hand side will contain an infinite number of terms, obtained from 
repeated reflection of R2 in the plates. We emphasize that we do not assume that these 
reflecting boundary conditions still hold in a system for which the plates have a prefer- 
ence for one of the monomer types over the other. We only assume that the Landau free 
energy of the general system can be written as a sum of two parts, one being the Lan- 
dau free energy of the system without weak interactions, and that reflecting boundary 
conditions may be assumed in the calculations concerning the latter. 

We introduce a coordinate system for which the z— axis is perpendicular to the film, 
and for which the origin is chosen such that the film is in between z = and z = d. 
The profile \|/(x,j,z) is only defined within the film; that is, for <= z <= d. We 
now construct the profile ^(x, y, z) , which is defined in the whole of space, by repeated 
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reflection of the film profile in the plates. It can be shown that under the assumptions 
enunciated above, the free energy density of a film with profile \\f{x,y,z) is the same as 
the free energy density of a bulk system with profile Since much is already 

known about the calculation of the free energy density of bulk systems, it is convenient 
to switch from the film to the auxiliary bulk system. The profile of the auxiliary bulk 
system is invariant under translation over a distance 2d into the z-direction, and mirror 
symmetric in the plane z = 0. Note that it follows from this that the profile is mirror 
symmetric in any plane described hy z = kd for integer k. The symmetries imply that 
we can express the composition profile as a Fourier series of the following general 
form: 



klZz 

^{x,y,z)=i:';:=iAk{x,y) cos— (4) 

For even values of k, the cosine appearing in equation ^ is symmetric with respect 
to the midplane of the film, while for odd values of k, it is antisymmetric. Since the 
plates are identical, only for even values of k does the corresponding composition wave 
couple to the plates. Consider for given vector q and phase (|) the following combination 
of plane waves: 



g-'<i> gi{-1xx-qyy + qzz) _j_ ^ii^ i{- qxX~ qyy + qzz) 

= 2 COS {qj^x + qyy + ^) cos{ q^z ) 

with 

= -J (6) 

We see that any real linear combination of profiles of the form given by equa- 
tions ^ and ^ has the special form given by equation Q, and thus possesses the 
symmetries of an auxiliary bulk system. In fact, the converse is also true: any bulk 
profile satisfying the above-mentioned symmetry conditions can be written as a linear 
combination of functions of the form given by equation 

= Y.A^^'^^^.^ix.y.z) (7) 

where the summation over q is restricted to some appropriate quadrant of Fourier 
space (see equation Q), and the amplitudes As are real. It is easier to express the 
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symmetry requirements in terms of the Fourier transform of the composition profile, 
which is defined by 

vj/. = y cPr ^(r) e'^-' (8) 

If ^(r) is periodic, then its Fourier transform consists of a discrete collection of 
5— peaks, each of which has an amplitude proportional to the volume of the system. 
Using equation (jSj), we see that in order for the profile to be real and the symmetry 
conditions to be satisfied, the following set of three conditions on the Fourier transform 
of the composition profile is sufficient (it is also necessary, but we will not prove this) 



^(-^) = ^*(^) (9) 

kit 

^(Q/j^Iz) = if not <lz = Y 

where k is an integer, the star denotes a complex conjugate, and q^^ = {q^.qy) 
denotes a vector parallel to the film. The first condition expresses the fact that the 
composition profile must be real-valued, while the second and the third express the 
symmetry properties. Note that it follows from the second condition that ^{q) is 
real for any vector which is perpendicular to the film, which means that sines are not 
allowed; only cosines are. For the set of ^—vectors which attain a non-zero amplitude 
the second condition implies that whenever [q^ , q,) is in the set, then (— , qz) is also 
in the set, and both have the same amplitude, but opposite phases. We stress that the 
symmetry conditions expressed by equation © are imposed on the profile ^ of the 
auxiliary bulk system, not on the profile \|/ of the film, which is the restriction of ^ to 
the film. The film profile is not constrained by symmetry requirements, though it is 
constrained by the requirements that it must integrate to zero, and that its value must 
lie between (— /) and (1 — /). 



3 Weak segregation regime 

We will study the film in the weak segregation regime, which is defined by the re- 
quirements that the A5— interaction strength % is close to its critical value %c, the 
A— monomer fraction / is close to 1/2, and the preference of the plates for one of 
the monomer types over the other is only slight. The precise meaning of the words 
"close" and "slight" will be given further on. We will also assume that the film thick- 
ness is close to a multiple hLq of the natural periodicity Lq = 2n/qo, where qo denotes 
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the position of the minimum of the second order vertex ["T| . In our calculations we will 
assume that the structure is periodic with period d/n'm the direction perpendicular to 
the film, by which we mean that if the film is divided into n slices of thickness d/n, 
then each slice is identical as far as the internal structure is concerned. This assump- 
tion of homogeneity is only justified if the film thickness d is much smaller than the 
bulk correlation length ^, which depends on the system parameters via ll23l 

where Rg is the radius of gyration of an unperturbed chain in the melt. If, on the 
other hand, the film thickness is much larger than ^, then in the middle of the film 
the structure will be the equilibrium bulk structure, which may be different from the 
structure near the plates. Nevertheless, for any given value of the film thickness d, no 
matter how large, the condition ^» d will be fulfilled once we are close enough to the 
critical point. Therefore, in the weak segregation regime our assumption concerning 
the homogeneity of the structure in the perpendicular direction does not impose any 
restriction on the film thickness, and the integer n may take on any value. For large n 
we are in the regime Lq << J << ^. Due to the first inequality, incommensurability 
effects are negligible. Due to the second inequality, the free energy contribution com- 
ing from the surface preference is spread evenly over the entire film, and since the film 
is thick while the surface preference is small, the effects of the surface preference will 
be negligible as well. We conclude that if n is large, neither the incommensurability, 
nor the surface preference plays an appreciable role, so that the structure inside the 
film will concide with the bulk structure. For this reason we need only to consider 
small values of n. We stress that for given film thickness d» Lq the reasoning given 
above only applies once we are close enough to the critical point, where ^» d. If, 
on the other hand, ^ < J, the influence of the surface preference penetrates only over 
a certain distance into the film, and the structure near the midplane may be different 
from the structure near the surfaces. It would be interesting to study how one structure 
slowly transforms into another; see also I.24J . where this problem was addressed in the 
limit d ^ °°. 



4 Free energy 

The next step is to write down the free energy of the auxiliary bulk system in terms 
of its composition profile m. Since we are in the weak segregation regime, the free 
energy may be expanded in powers of the composition profile. It can be shown that in 
the generic case, the terms till the fourth order have the same order of magnitude, while 
the higher order terms are negligible. This leads to the following approximation to the 
free energy of the auxiliary bulk system (we did not yet include a term accounting for 
the preference of the plates for one of the monomer types): 
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The interaction between A-monomers and B-monomers has been taken into ac- 
count in the usual way by means of the Flory-Huggins %— parameter. If the composi- 
tion profile is periodic, the Fourier transform is only non-zero on a discrete lattice, 
called the reciprocal lattice. In the weak segregation regime, the spacing of this lattice 
is given by the position of the minimum of the second order vertex function F^^). 
The lattice vectors having length ^o, forming what is called the first harmonics sphere, 
attain a much larger amplitude than the vectors of the higher harmonics spheres. For 
this reason the latter are usually neglected. In this so-called first harmonics approxima- 
tion all vectors appearing in the summations of equation (fTTT) have the same length qo. 
Strictly speaking, this can only be arranged if the film thickness is a multiple of half 
the natural periodicity (and even in that case not all structures can be made to fit; see 
the next section). For other values of the film thickness, the Fourier vectors which are 
perpendicular to the film (which will always be present if the plates have a preference 
for one of the monomer types over the other) will have to be compressed or extended, 
giving them a length different from qQ. However, since we will only consider values 
for the film thickness which are close to a multiple of the natural periodicity, the length 
of these vectors will still be close to qo. The fact that all vectors appearing in the sum- 
mations of equation (fTTT) have more or less the same length can be used to simplify 
this equation. The third and fourth order vertices F^^) and F^^^ depend only weakly on 
the angles between their arguments, and may be approximated by angle-independent 
constants, as follows: 

The constants ^ and X will depend on the composition /. For the second order ver- 
tex this approximation would be too crude, because although the length of all Fourier 
vectors will be close to qo in the regime under study, the relative change in the value of 
the second order vertex is large if we change the length of the ^—vector even slightly, 
because F*^^^ (q) is close to zero for q = qo (where q denotes the length of the vector 
q). Nevertheless, it is not necessary to take into account the full ^—dependence of the 
second order vertex; it will be sufficient to make a quadratic approximation around its 
minimum q = qo, as follows: 

F(2)(^)-2x ^ c{q-qo)^-X (13) 
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The parameter x represents the distance to the spinodal, and it is a function both 
of X, and of /. The free energy density of the film can be obtained from that of the 
auxiliary bulk system by adding a term Opiates accounting for the interaction of the 
monomers with the plates (see also equation |2l). Note that we are talking about the 
weak interaction related to the preference of the plates for one of the monomer types 
over the other, not the excluded volume interaction. Under the assumption that this 
weak part of the monomer-plate interaction is pair- wise additive and short-ranged (the 
range should be much smaller than the radius of gyration of the blocks), the corre- 
sponding contribution to the free energy coming from one of the plates will be linear 
in the A— monomer fraction near this plate, averaged over its 2-dimensional surface. 
We thus obtain for some constant Hq 



^plates = -/lo(^(^,y,0)),,, - /lo(^(x,y,J)).,, (14) 

We will not attempt to express ho in terms of more basic quantities, but take it 
as a phenomenological parameter. Adding the various terms together we obtain the 
following expression for the free energy: 



XF 1 ^ (2) ^ . £ \-i 

^ ^ Q {Qi} 



+ i E -^21^22^23-^24 - ^ E*^2 (15) 
{Qi} ^ 2 



where Qi for / = 1,2,3,4 represents a lattice vector, and the summations are re- 
stricted to sets of lattice vectors whose sum equals zero. The star in the last summation 
indicates that it is restricted to lattice vectors which are perpendicular to the film, and 
whose length satisfies equation ® for even values of k. The variables xq are rescaled 
amplitudes defined by (see equation O for the (implicit) definition of Aq) 



A,^^ ,16, 
The constant h is defined in terms of ho (see equation (IT?!) ') by 



h = ^ (17) 

where Lq = 2n/ qo is the natural periodicity. The parameter his a measure for the 
preference of the plates for the monomers of type A (if h > 0) or B (if h < 0). Note 
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that in the last term of equation (fTTT) the division by d is needed in order to arrive at the 
free energy density, since the volume of the film is proportional to d. The parameter 
£ appearing in equation dTSt has been defined in such a way that it is positive if the 
A— monomers are in the minority. It is given by (see equation (fT2ll ') 



The parameter £ vanishes for / = 0.5, and close to the critical point its dependence 
on the composition / and the chain length N may be approximated by 



where ^ is a constant. 



^ if-h (19) 



5 Microstructures in thin films 

We have remarked in section Q that if /z = 0, then in the weak segregation regime the 
free energy density of a film equals the free energy density of the bulk system whose 
structure is obtained by repeated reflection of the film in the interfaces. Suppose that 
for a given film thickness there exists a film structure which would, upon repeated re- 
flection in its interfaces, lead to the equilibrium bulk structure. It follows that this film 
structure would minimize the film free energy. Such film structure has to be a slice of 
the equilibrium bulk structure, but the opposite is not true: not every slice of the equi- 
librium bulk structure whose thickness coincides with the film thickness would lead 
back to the equilibrium bulk structure upon repeated reflection in its interfaces. For a 
given film thickness d such slices may not even exist, depending on the periodicity Lq 
and the symmetry of the equilibrium bulk structure. There is no general rule involv- 
ing the periodicity; for instance, the condition d = kLo for some positive integer k, is 
neither necessary nor sufficient. The fact that this condition is not sufficient becomes 
clear if we consider the gyroid structure. Its weak segregation Fourier transform con- 
tains twelve plane waves pointing into various directions, so that in any direction the 
gyroid structure is either non-periodic (but non-constant), or it has a period which is 
larger than twice the natural periodicity. It follows that there is no symmetric slice 
whose thickness equals the natural periodicity. The fact that the above-mentioned con- 
dition is not necessary becomes clear if we consider the lamellar structure. For a given 
periodicity, any slice which is cut perpendicular to the lamellae would lead back to 
the same lamellar structure upon repeated reflection in the plates, irrespective of the 
thickness of the slice. 

In order to avoid cumbersome formulations, we will say that a slice of a bulk struc- 
ture is 'symmetric' if it leads back to the same bulk structure upon repeated reflection 
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in its interfaces. Using this terminology, we can summarize the above by saying that 
if the equilibrium bulk structure allows for a symmetric slice having thickness d, then 
the equilibrium film structure for = is given by the structure inside the slice. In this 
form the principle is not very powerful, not only because its applicablity is restricted 
to the case h = 0, but also because symmetric slices exist only for a discrete set of 
film thicknesses. However, it does suggest a way to search for the equilibrium film 
structure in the more general situation. The idea is as follows. First we choose a bulk 
structure, and we try to find a symmetric slice whose thickness is close to the film 
thickness. We compress or expand the bulk structure in the direction perpendicular to 
the slice until the thickness of the slice coincides with that of the film. We may also 
want to compress the structure in some direction(s) parallel to the film in order to reset 
the lengths of the oblique Fourier vectors back to their natural value qo. Finally we 
calculate the free energy of the resulting deformed structure, adding a term to account 
for the surface preference. We repeat this procedure for various bulk structures, and 
compare their free energies. 

We have to search for symmetric slices of small integer thickness n (from now on 
we will measure lengths in units of the natural periodicity Lq = In/qo). This search is 
facilitated by considering the bulk structures in Fourier space. The Fourier transform 
of a periodic structure is only non-zero on a discrete set of Fourier vectors, which is 
called the reciprocal lattice. In a previous section we have given a symmetry condition 
on the reciprocal lattice in order for the associated structure to have a symmetric slice 
whose lower interface coincides with the plane z = 0. This symmetry condition is 
that if {qji^q^) is in the reciprocal lattice, then {—q//,qz), {q//, —Qz)-, and {—q//,—qz) 
must also be in the lattice. In addition, the z— component of any wave vector must be a 
multiple of n/d. Since the film thickness is given hy d = uLq, where Lq = In/q^, the 
z— components must be multiples of qo / In. For instance, if the slice has unit thickness, 
there only three allowed values for the angle a between the vector and the film, namely 
a = 0, a = 7r/6, and a = n/l. 

We will search for symmetric slices in the lamellar structure, the hexagonal struc- 
ture, the bcc structure, and the fee structure, which comes down to finding orientations 
of their respective reciprocal lattices satisfying the above-defined conditions. If a sym- 
metric orientation does not have perpendicular vectors, we will have to add them by 
hand, because such vectors will always have a non-zero amplitude if the plates prefer 
one of the monomer types over the other. 

We first consider the lamellar structure. In the first harmonics approximation, its 
Fourier transform consists of one pair of opposite vectors. There are only two ways 
of satisfying the symmetry conditions formulated above: the vectors should be either 
perpendicular, or parallel to the film; see Figure ([l]). Note that in the latter case we 
have added by hand a set of perpendicular vectors. In real space, the first Fourier 
transform corresponds to lamellae which are oriented parallel to the plates, while the 
second Fourier transform corresponds to lamellae which are oriented perpendicular to 
the plates. The perpendicular Fourier vectors make it possible that the lamellae have a 
different thickness near the plates as compared to near the midplane. Figure Q shows 



11 



a cross section of the parallel structure, while Figure Q shows a cross section of the 
perpendicular structure, both for the case n = \. Figure ® shows a cross section of 
the perpendicular structure for the case n = 2. Note that there is a thickening of the 
minority phase near the midplane, which is due to the natural periodicity of the block 
copolymer. 

Next we consider the structure which consists of hexagonally packed parallel cylin- 
ders. In the first harmonics approximation, its Fourier transform is a set of six coplanar 
vectors of length qo, making angles of sixty degrees with each other. There are three 
ways of orienting the vectors with respect to the film such that the symmetry conditions 
are satisfied, which are illustrated in Figure Q- The first possibility is to orient all six 
vectors parallel to the film, in which case we need an additional set of perpendicular 
vectors. In real space this corresponds to hexagonally packed cylinders oriented per- 
pendicular to the plates. Due to the perpendicular Fourier vectors, the cylinders will 
have a different thickness near the plates as compared to near the midplane (see Figure 
The second possibility is to orient two of the six vectors vectors perpendicular to 
the film. Since the z— components of the four oblique vectors are ±qo/2, this structure 
will fit without distortion in any film whose thickness is a multiple of the natural peri- 
odicity Lq. In real space it corresponds to the cylinders lying parallel to the film, and 
a cross section for the case n = 1 is given by Figure dTj. The third way to orient the 
structure is obtained if we rotate the previous (parallel) structure over an angle of 30 
degrees, so that two Fourier vectors become parallel to the film. We will have to add 
a set of perpendicular vectors to account for the plate interaction. The z— components 
of the oblique vectors are equal to ^/3/2 in units of qo, and since this is an irrational 
number, the structure will not fit in any film of integer thickness. If n is large enough, 
it can be made to fit by just a slight distortion, but the corresponding free energy will 
always be higher than that of the perpendicular orientation, since both structures have 
the same terms in the free energy. It follows that this way of orienting the cylinders is 
never stable, which is in accordance with experimental findings [ilOil . 

Now we consider the bcc structure. Apart from trivial rotations around a perpen- 
dicular axis, there are only two possiblities to orient the Fourier vectors such that the 
symmetry requirements are satisfied (see Figure @). The first possibility is to orient 
one Fourier vector perpendicular to the film. All oblique vectors have a z— component 
which is equal to qQ/2, and the structure fits without distortion in any film of integer 
thickness. Figure ® shows what this structure looks like in real space in a film of unit 
thickness. If the film thickness is slightly increased from integer value, the oblique 
vectors have to be turned slightly towards the horizontal orientation in order for their 
z— component to remain commensurate with the film thickness. At the same time, the 
angle between two oblique vectors whose sum equals one of the horizontal vectors 
must be increased. Since we ignore the angle dependence of the vertices, these dis- 
tortions do not complicate the calculations. In the second way to orient the Fourier 
vectors, none of them is perpendicular to the film, and perpendicular vectors have to 
be added in order to account for the plate interaction. The z— component of the oblique 
vectors equals q^j \f2, and since \/2 is irrational, the structure will not fit without dis- 
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tortion in any film of integer thickness. This means that for small values of n this 
structure would have a high free energy and may be ignored, but for large values of n 
it would have to be taken into account, since the degree of distortion needed to make 
it fit would be small. In fact, it can be argued that it will become stable in thick films 
if the plates have a preference for the majority component. Nevertheless, for reasons 
expounded before, we restrict ourselves to small values of n, so that we may neglect 
the second orientation of the bcc structure. 

There is also a film structure which is derived from the bulk fee structure. We 
took it into account in our calculations, but since it is nowhere stable, we will refrain 
from discussing it any further. Finally we considered the gyroid structure. Due to its 
large number of Fourier vectors pointing into various directions, it does not fit inside 
thin films (i.e., n small), and can only be embedded in thick films. Since we restrict 
ourselves to small values of n, we do not need to consider the gyroid structure. 

6 Rescaled parameters 

Using the results of the previous section, we can write down expressions for the free 
energy of the various structures. These expressions will not only depend on the sym- 
metry of the structure and the amplitudes of the Fourier vectors, but also on the degree 
of distortion (compression or stretching), which can be quantified by the dimensionless 
parameter D defined by 

d 

D = — (20) 
nLo 

where the integer n is fixed by the requirement that nLo is as close as possible to d. 
The structure is compressed if D < 1, and stretched if D > 1. Since we consider films 
whose thickness is close to a multiple of the natural periodicity, we have D ^ I. We 
are now in a position to write down the expressions for the free energy of the various 
structures in terms of the Fourier amplitudes defined by Figures ([T]), Q and Q. The 
result is 
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A- ~9 9 ^ d ^ A. 79 

TjFsqr = -XM^ - XV^ + -a* + -V^ + mV - 

K 4 4 nL) 

— Fhex = - 2XV^ - 2£MV'^ + -u' + -v^ + 2mV 

y // 4 2 nD 

^ „ ~9 ^9 IzL 15/1 ^99 

TjFhex, = -XM^ - 3xv^ - 2£v^ + -a^ + — + 3mV - 

y 1 4 4 nD 

^ ~9 9^9^ 9^ 9I4.I4. 

— Ffe^^ -XM - XV - 4xw - 4euw - 4evw + -u + -v 

+ Qvv"^ + M^v^ + 4m2w^ + 4v^w^ + 4mvw^ ^ (21) 

nD 

The parameter x has been defined in equation (fT^). and is given by 

-X = r(2)(<5ro)-2% (22) 

while X is defined by 



- X 



r(2)(^o/^)-2x (23) 



By making use of the fact that D is close to unity, we can simplify the free energy 
expressions of equation (OTT) . reducing the number of parameters from four (namely, /, 
X, h, and D), to three. First we relate x to x by combining equation (fOl with equation 



-X = cql(^^^ -X ~ cql{D-\f-x (24) 

where we have only kept the leading order term in (1 — D), which is justified by 
our assumption that D is close to unity. In the same approximation, the factor D may 
be omitted from the linear term describing the interaction between the A-monomers 
and the plates; that is, in equation (OTT) we may make the substitution 



2hu 2hu 
nD n 
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As a first step in the reduction of the number of parameters, we replace the original 
set by the equivalent set (£,x,/z,Z)). Next we rescale the parameters x, h, 

and D — 1 by some appropriate power of £, and arrive thus at the new parameters A, T, 
and//: 



^ |D-1| X h 

A=^qoL^\ T= //= (26) 

|£| 

We also rescale the amplitudes 

u ^ £u v^ev w ^ £w (27) 

and the free energy 



Substitution of equations (l26t . (l27t . and (l28t into equation (|2T1) leads to 



(28) 
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-T)u^ 
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+ 4" - 
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- Tv^ + 


1 4 22,14 2Hu 

-U^ + U V + -V 

4 An 


Fhex 1 J 


= (A^- 


-T)u^ 


- 2Tv^ - 


4 2 n 


Fhex ^ 


= (A^ 


-T)u^ 


- 3Tv^ - 


_T I4 15 4 2//m 

- 2v^ + -M^ + 3m^v^ + — 

4 An 


Fbcc 


= (A^ 


-T)u^ 


-Tv^- 


ATw^ -Auw^ -Avw^ + + -v^ + 

4 4 



+ 9^"^ + M^v^ + 4m2w^ + Av^w^ + 4mvw^ - 

n 

(29) 

From the fact that the (small) parameters x, £, h, and D—\ appear only in the com- 
binations r, //, and A, it can be concluded that the relevant values for these parameters 
satisfy the scaling 



x~£^ (D-l)~£ h^z^ (30) 
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We note that although the original parameters (/ — 0.5), (% — Xc), {D — I), and h 
all have to be small in order for the theory to be valid, expression (1^^ is valid for 
any values of the parameters A, T, and H. This can be understood by realizing that 
we went from four independent parameters to three. There must be a fourth indepen- 
dent parameter in addition to A, T, and H, say X, on which the free energy depends 
in a trivial manner, so that the ratio of any two free energies does not depend on X. 
This parameter X is a measure for the distance to the critical point; therefore, it must 
be small in the weak segregation regime. The remaining three parameters describe the 
position on the three-dimensional hyper-surface defined by the condition that X is con- 
stant. Varying these parameters between zero and infinity only changes the position on 
this hypersurface, without bringing us further away from the critical point. It follows 
that in the phase diagrams to be presented in the next section, where A is placed along 
the horizontal axis and T along the vertical axis, any region in the phase diagram is 
meaningful in the weak segregation regime, not just the region near the origin. 

7 Results and Discussion 

We will present our results in the form of phase diagrams. Along the horizontal axis we 
put A ~ |D — 1 1/|£|, and along the vertical axis we put T = %/£?-. Since there is a third 
parameter, namely H = h/£?, we need to calculate a series of such phase diagrams, 
for increasing values of H. From the definitions of the parameters A, T , and H we see 
that the phase diagram does not change if we flip simultaneously the signs of £ and 
h. This is due to the fact that if we exchange A and B in the copolymer molecules, 
and at the same time reverse the plate preference, the equilibrium structure is simply 
inverted according to \j/ ^ — \|/, keeping the same symmetry and amplitude. Therefore, 
the phase diagrams for // > describe both the situation where the A— monomers 
are in the minority (£ > 0) and the plates prefer the A— monomers {h > 0), and the 
situation where the 5— monomers are in the minority (£ < 0) and the plates prefer the 
5— monomers (h < 0). In order to avoid cumbersome formulations, in our discussions 
we will always assume that the A— monomers are in the minority. 

The results are shown in Figures (fTTTl to (l20b . We will start with the simplest case 
H = (Figure [TSt. where the plates have no preference for either monomer type. If 
the film thickness is an integer, corresponding to A = 0, the phase behavior is identical 
to that of a bulk system, because the weak segregation bulk equilibrium structures all 
fit without distortion inside a film of unit thickness. We see that on increasing the 
parameter T, which is a measure for the A5— repulsion, the order in which the phases 
appear is given by homogeneous - bcc - hexagonal - lamellar, which is in accordance 
with well-known results |3|. If the film thickness is not an integer, so that A > 0, the 
lamellae and the cylinders are oriented perpendicular to the film, and the free energy 
of the corresponding structures is independent of the film thickness. This explains why 
the phase boundary separating them is a straight horizontal line. The free energy of 
the three-dimensional bcc structure, on the other hand, increases if the film thickness 
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deviates from being an integer, because it will be compressed or stretched. This is the 
reason why the bcc structure disappears upon moving to the right in the phase diagram. 

Next consider the phase diagrams for positive values of H. We see that close to 
the line A = 0, the lamellae and the cylinders have a parallel orientation. The reason 
is that in this way the film can profit optimally from the preference of the plates for 
the monomers of type A (remember that without affecting the generality we will as- 
sume that the A— monomers are in the minority). It is for this reason that the border 
between the parallel and the perpendicular phases shifts to the right if H is increased. 
Although the perpendicular hexagonal structure also has a set of perpendicular Fourier 
vectors which couple to the plates (see Figure 0), these vectors couple to the other 
Fourier vectors only via the (positive) fourth order term, while the perpendicular vec- 
tors of the parallel hexagonal structure couple favorably to the other Fourier vectors 
via the (negative) third order term. The stability of the parallel lamellar structure over 
the perpendicular lamellar structure for small values of A can also be explained by 
the positive fourth order contribution to the free energy in the latter coming from the 
coupling between the parallel and the perpendicular Fourier vectors. The stronger the 
plate preference H, the more the parallel structures are favored over the perpendicular 
structures. More generally, the shifts in the phase boundaries for increasing value of 
H can be accounted for by differences in the amplitude of the perpendicular Fourier 
mode. The structure with the largest amplitude will profit most from an increase in the 
value of H (note that the amplitude may be negative). 

Next we consider the phase diagram for H < 0. The most striking feature is the 
total absence of the parallel hexagonal phase. This is due to the fact that if // < 0, 
the third order term and the first order term of the free energy have opposite signs, 
and compete with each other (see equation 1^^. In case the minority component forms 
the cylinders while the majority component forms the matrix, the third order term 
is negative, but there will be an excess of the minority component near the plates, 
which is unfavorable for H <0. If, on the other hand, the structure is inverted so that 
the majority component forms the cylinders while the minority component forms the 
matrix, the former will be in excess near the plates, which is favorable, but the price 
to be payed is a positive third order contribution. Due to this frustration the parallel 
phase is nowhere stable for H < 0. One might think that it should be possible to cut 
a slice out of the hexagonal bulk structure such that the cylinders are in the middle 
of the film, so that the material making up the matrix is in excess near the plates, 
preventing the frustration between the first and the third order terms of the free energy 
expansion. However, it appears that any slice of integer thickness having the cylinders 
in the middle violates the symmetry conditions; see, for instance, the failed attempt 
illustrated by Figure (flOb . The disappearance of the bcc structure for negative values 
of H can be accounted for by the same mechanism. It should however be noted that 
for large values of n, the alternative way of cutting a slice out of the bcc structure, as 
illustrated by the second picture of Figure (), would lead to a lower free energy, so that 
the corresponding bcc2 structure might remain stable even for negative values of H. 
Finally we remark that for any point (A, %) in the phase diagram, the parallel lamellar 
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structure becomes stable once the absolute value of H is large enough. 

Quantitative comparison with experimental work or earlier theoretical work is not 
well possible, but some of the qualitative aspects of our phase diagrams are in accor- 
dance with earlier findings. In |20| it was remarked that if the surface potential H 
increases, the parallel structures are favored over the perpendicular structures, which 
is confirmed by our series of phase diagrams. The reason can be found in the fact that 
in the parallel structures, the excess of the minority component at the plates is larger 
than in the perpendicular ones, while the minority component is favored by the plates 
for positive H. 

In [ 19 1 it was observed that the lamellar structure gains stability over the hexagonal 
structure if the surface potential increases, which is also apparent from our series of 
phase diagrams. 

In [17 1 a symmetric diblock copolymer trapped between two parallel plates was 
investigated using both self-consistent field theory, and Monte Carlo simulations. It 
was found that in thick films, the lamellae are always oriented parallel to the plates, 
while for thin films there are transitions between parallel and perpendicular orienta- 
tions. Strictly speaking, our phase diagrams are not valid for symmetric copolymers, 
since we have rescaled our parameters with powers of £. Of course, the analysis would 
be easier for symmetric copolymers, and one could represent their phase behavior by a 
single diagram in which x//i^/^ is plotted versus \D— \ \/h^^^ . However, it is possible 
to understand the occurrence of periodic transitions between parallel and perpendicu- 
lar orientations of the lamellar structure on the basis of our phase diagrams. Suppose 
that for a given copolymer we slowly increase the film thickness, starting from unity. 
What will be the corresponding path in our phase diagrams? For unit film thickness 
we are somewhere on the vertical axis, and we move to the right if the film thickness 
is increased. However, if we increase the film thickness further, at a certain moment 
the number of periods fitting inside the film will increase jumpwise (either two 1 .5, or 
to 2), with as a result that suddenly the structure is compressed rather than stretched. 
Increasing the film thickness further wil now lead to a decrease of the parameter D (see 
the previous section), so that we start moving to the left in the phase diagram, until we 
hit the vertical axis once again. It follows that if we keep increasing the film thickness, 
we will periodically move to the right and to the left. The turning point at the left will 
always lie on the vertical axis, but the turning point at the right will move to the left 
for thicker films, because the degree of frustration in thick films cannot be very large. 
Therefore, for thick films the amplitude of the periodic motion in the phase diagram 
is so small that the phase boundary between the parallel phase and the perpendicu- 
lar phase is never reached, so that the system stays permanently in the parallel phase. 
For thinner films, however, there might be periodic transitions between the two phases 
upon increasing the film thickness. We conclude that our results are in agreement with 
the results of reference lITTll . 
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Figure 1: Fourier vectors and amplitudes: lamellar structures 
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Figure 2: Fourier vectors and amplitudes: cylindrical structures 
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BCC 1 BCC 2 

Figure 3: Fourier vectors and amplitudes: bcc structure 




Figure 5: Lamellar perpendicular, side view, n=l 





Figure 7: Hexagonal parallel, side view, n=l 




Figure 8: Hexagonal perpendicular, n=l 
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Figure 9: BCC top view, n=l 




Figure 10: Slice through hexagonal phase violating symmetry conditions 
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Figure 11: Phase diagram for H = —0.80. 

11= lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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Figure 12: Phase diagram for // = —0.40. 

11= lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 



H = -0.2 




Figure 13: Phase diagram for // = —0.20. 

11= lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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Figure 14: Phase diagram for // = —0.05. 

11= lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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Figure 15: Phase diagram for // = 0. 

= homogeneous, 1 = lamellar, 2 = hexagonal, 3 = bcc 
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Figure 16: Phase diagram for // = 0.05. 

11 = lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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H = 0.1 




Figure 17: Phase diagram for // = 0. 1 . 

11 = lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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H = 0.2 




Figure 18: Phase diagram for // = 0.20. 

11 = lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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Figure 19: Phase diagram for // = 0.40. 

11 = lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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H = 0.8 




Figure 20: Phase diagram for H = 0.80. 

11 = lamellar parallel, 12 = lamellar perpendicular, 

21 = hexagonal parallel, 22 = hexagonal perpendicular, 3 = bcc 
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